function y = phi0phi0(x)
    y = phi0(x) .* phi0(x);
end